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Abstract 

We have considered nonuniversal gaugino mass models of supergravity, arising from 
a mixture of two superfield contributions to the gauge kinetic term, belonging to a 
singlet and a nonsinglet representation of the GUT group. In particular we analyse two 
models, where the contributing superfields belong to the singlet and the 75-dimensional, 
and the singlet and the 200-dimensional representations of SU(5). The resulting lightest 
superparticle is a mixed bino-higgsino state in the first case and a mixed bino-wino- 
higgsino state in the second. In both cases one obtains cosmologically compatible dark 
matter relic density over broad regions of the parameter space. We predict promising 
signals in direct dark matter detection experiments as well as in indirect detection 
experiments via high energy neutrinos coming from their pair-annihilation in the Sun. 
Besides, we find interesting 7-ray signal rates that will be probed in the Fermi Gamma- 
ray Space Telescope. We also expect promising collider signals at LHC in both cases. 
PACS No: 04. 65. -he, 13.40Em, 14.60Ef, 13.85.-t, 14.80.Ly 

1 Introduction 



A leading candidate for the cold dark matter of the universe is the lightest supersymmetric 
particle (LSP) of the minimal supersymmetric standard model (MSSM) [1]. Astrophysical 
constraints require the LSP to be colorless and chargeless, while direct dark matter (DM) 
search experiments strongly disfavour a sneutrino LSP. Thus the favoured candidate for LSP 
in the MSSM is the lightest neutralino 

X = x'i = ciB + C2W' + o,Hl + C4^^}. (1) 
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The neutralino mass matrix is given by 
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where Mi,M2 and /i are the bino,wino and higgsino mass parameters; and tan/? represents 
the ratio of two Higgs vacuum expectation values. The LSP is the hghtest eigenstate of the 
mass matrix. It can be a dominantly bino,wino or higgsino state or else a large admixture 
of these interaction eigenstates. 

Much of the MSSM phenomenology so far has been done in the context of minimal 
supergravity (mSUGRA) [2] model, because of its simplicity and economy of parameters. 
This model predicts a dominantly bino LSP over the bulk of the parameter space. The bino 
does not couple to gauge bosons since it does not carry any gauge charge. So the main 
annihilation mechanism of the bino dark matter is via t-channel slepton exchange, XX ~^ 
However the large slepton mass limits from LEP [3], particularly in the mSUGRA 
model, makes this annihilation mechanism inefficient. This leads to a gross overabundance of 
the DM relic density over the bulk of the mSUGRA parameter space. There are only narrow 
strips of parameter space, like the stau coannihilation (xf ^ t'j) and resonant annihilation 
{XX ~^ bb,ti,T~^T~) regions giving cosmologically compatible DM relic density [4] which 
require however some stringent mass correlations between the annihilating LSPs and the 
intermediate particle in the s-channel resonance or between a coannihilating sparticle and 
the LSP. 

In contrast, the higgsino and wino carry weak isospin I = \ and 1 respectively. So a 
higgsino and wino-dominated LSP can pair-annihilate efficiently via their couplings to the 
W/Z gauge bosons, leading to an underabundance of DM relic density for sub-TeV LSP 
masses. One can get cosmologically compatible relic density only for relatively large LSP 
masses of Mj^ ~ 1 TeV and ^ 2 TeV, which make these models inaccessible to LHC. 
They can only be probed at a multi-TeV linear collider like CLIC [5,6]. 

It is evident from the above discussion that a mixed bino-higgsino or bino-wino LSP 
is expected to give cosmologically compatible DM relic density for sub-TeV LSP masses, 
which can be probed at the LHC. This has been described as the 'well-tempered' scenario 
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in Ref. [7]. The DM relic density and the detection phenomenology of this scenario have 
been studied in a model independent way in Ref. [8]. It should be mentioned here that a 
very important example of the mixed bino-higgsino LSP occurs in the so-called hyperbolic 
branch/focus point (HB/FP) [9-12] region of mSUGRA model. This is a narrow strip at 
the edge of the mSUGRA parameter space, corresponding to very large (multi-TeV) scalar 
masses. 

In this work, we investigate a class of nonuniversal gaugino mass models, where the mixed 
neutralino LSP can be realised in a simple and economical way. In particular we construct 
two models — the first giving a mixed bino-higgsino LSP and the second giving a mixed 
bino-wino LSP. The latter model also gives a triply mixed bino-wino-higgsino LSP over a 
significant part of the parameter space. For each of these two models we find cosmologically 
compatible DM relic density over broad bands passing through the middle of the parameter 
space in the mo — mi plane. We also investigate the expected signals in both cases for direct 
and indirect DM detection experiments. In each case we find promising signals for future 
direct detection experiments. Likewise we find promising indirect detection signals in each 
case for the IceCube experiment, in the form of high energy neutrinos coming from DM 
pair-annihilation inside the Sun. We also estimate the corresponding line and continuum 
7 ray signals coming from their pair-annihilation in the galactic core. We conclude with a 
brief discussion of the expected collider signatures of these models at LHC 



2 Nonuniversal Gaugino Mass Models for mixed neu- 
tralino DM: 

SUGRA models with nonuniversal gaugino masses at the GUT scale have been studied in 
many earlier works [13-21]. We shall only summarize the main results here, focusing on the 
simplest and most predictive GUT group, SU{5). Here the gauge kinetic function that is 
related to the GUT scale gaugino masses, arises from the vacuum expectation value of the 
F-term of a chiral superfield which is responsible for SUSY breaking. This results in a 
dimension five term in the lagrangian, 

L D ±£^X^X.^ (3) 

J^Planck 
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where Ai,2,3 are the f/(l), SU{2) and SU{3) gaugino fields — bino (B), wino (W) and gluino 
(g). Since the gauginos belong to the adjoint representation of SU{5), $ and F$ can belong 
to any of the irreducible representations appearing in their symmetric product, i.e. 

(24 X 24),^^^ = 1 + 24 + 75 + 200. (4) 

The minimal SUGRA model assumes $ to be a singlet, implying equal gaugino masses at the 
GUT scale. On the other hand if $ belongs to one of the nonsinglet representations of SU (5), 
then these gaugino masses are unequal but related to one another via the representation 
invariants. Thus the three gaugino masses at the GUT scale in a given representation n are 
determined in terms of a single SUSY breaking mass parameter m^^^ by 

^1^2,3 = ^^2,3''^l/2; (5) 

where 

Cl,^, = (1, 1, 1), C5,3 = (-1, -3, 2), C5,3 = (-5, 3, 1), C5°3 = (10, 2, 1). (6) 

The nonuniversal gaugino mass models are known to be consistent with the observed 
universality of the gauge couplings at the GUT scale [13,15], with 1/25). Since 

the gaugino masses evolve like the corresponding gauge couplings at one-loop level of the 
renormalisation group equations (RGE), the three gaugino masses at the electroweak (EW) 
scale are proportional to the corresponding gauge couplings, i.e. 

Ml = {ai/aG)M^ ~ (25/60)Cfm^/2 
M2 = {a2/aG)M^ ~ (25/30)C^m^/2 

M3 = {a3/aG)M^ :^{25/9)C^m^y,. (7) 

For simplicity we shall assume a universal SUSY breaking scalar mass mo at the GUT scale. 
Then the scalar masses at the electroweak scale are given by the renormalisation group 
evolution formulae [22,23]. A very important SUSY breaking mass parameter at this scale 
is Tuff^, which appears in the EW symmetry breaking condition, 

o o mir — mir tan^ n 
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The last equality holds for the tan/3 > 5 region, which is favoured by the Higgs mass limit 
from LEP2 [24]. Expressing mfj^ at the right in terms of the GUT scale parameters via the 
one-loop RGE gives [23] 



at a representative value of tan (3 = 10, neglecting the contribution from the trilinear coupling 
term at the GUT scale. Moreover, the coefficients vary rather mildly over the moderate tan f3 
region. 

Although we shall use exact numerical solutions to the two-loop RGE in our analy- 
sis, one important result is worth noting from the simple formulae of Eqs. (5,6,7,9). For 
the singlet and 24-plate representations of the SUSY breaking chiral superfield $, we have 
\Mi\ < IM2I, This corresponds to a bino LSP, resulting in overabundance of DM relic 
density over most of the parameter space. On the other hand, for the 75-plate and 200-plate 
representations we get < \Mi\, IM2I. This corresponds to a higgsino LSP, resulting in 
underabundance of DM relic density [17]. 

In general the gauge kinetic function may contain several chiral superfields, belonging to 
different representations of SU (5), which gives us the freedom to vary the relative magnitudes 
of Mi,M2 and /i continuously. In particular we shall consider two models, where the gauge 
kinetic function contains a mixture of singlet plus 75-plet superfields (1 -|- 75) in the first 
case and a singlet plus 200-plet superfields (1-1-200) in the second. By adjusting the relative 
contribution from the two superfields we shall obtain a mixed bino-higgsino LSP in the first 
case and a mixed bino-wino LSP in the second, each giving favourable DM relic density over 
large regions of the parameter space. It should be mentioned here that mixtures of singlet and 
nonsinglet superfield contributions were first discussed in Ref. [13]. Moreover, the (1 -|- 24), 
(1-1-75) and (1-1-200) models were considered in the context of i) direct detection of neutralino 
DM in Ref. [14], and ii) a bino DM in Ref [18], where in the latter the relative contribution 
from the two superfields were adjusted to recover the bulk annihilation region of the right 
DM relic density in each case. Nonuniversal gaugino mass models based on such mixed 
representations have also been considered more recently to investigate the Higgs production 
phenomenology via SUSY cascade decay [20]. To the best of our knowledge however, this 
is the first investigation of such nonuniversal gaugino mass models, to obtain the right relic 



-O.lm^ + 2.1Mf - Q.22Mf - 0.006Mf ^ + 0.006Mf + 
Q.l^M^M^ + 0.03Mf Mg^, 



(9) 
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density via mixed neutralino DM. As in [18], we shall vary the relative contribution of the 
two superfields via a single parameter, i.e 



^ - Q - a75)"^i/2 and ml% = "751111/2, (10) 



1/2 — l-L ~ LJ;75;/;ti/2 ciiiu iii]^/ 



for the (1+75) model and 



1/2 = (1 - «20o)"^i/2 and m^^^ = a2oomi/2, (11) 

for the (1+200) model. Then one sees from Eqs. (5,6,7,9) that for the (1+75) model, 

= 0.50 =^ \Mi\ ~ 1^1 ~ mi/2 < IM2I (12) 

for mo > mi/2- This leads to a mixed bino-higgsino LSP. One also sees from Eq. (9) 
that decreasing tuq leads to increase of and hence decrease of higgsino fraction of LSP. 
Moreover one sees from eqns (5,6,7) that for (1+200) model 

0200 = 0.1 ^ |Mi| ^ IM2I, (13) 

leading to a mixed bino-wino LSP. Again one sees from equation (9) that increasing mo leads 
to decrease of resulting in a triply mixed bino-wino-higgsino LSP in the large mo region. 
We shall present the results for the (1 + 75) and (1 + 200) models in the following sections. 
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0.55 
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448 


0.45 


1073 


584 
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Table 1: Variation of mass parameters with 075 (for the (1+75) model) and 0200 (for the 
(1+200) model) for mo = mi/2 = 500 GeV, Aq = 0,tan/9 = 10, and sign(/i) = +ve. 



We have probed how the nature of the LSP varies with the choice of mixing parameters 
in Table 1. We show Mi and fi for a few values of 075 for the (1+75) model for a given 
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set of other parameters. Similarly we show Mi, M2 and ji for a few values of 0200 for the 
(1+200) model. One can infer the dominant component of the LSP from the relative values 
of these masses. In the (1+75) model, the LSP is dominantly bino for 075 < 0.45, leading to 
a generic overabundance of DM. On the other hand it is dominantly higgsino for 075 > 0.55, 
leading to a generic underabundance of DM. In the (1+200) model, the LSP is dominantly 
wino for 0200 > 0.15, becoming dominantly higgsino for 0200 > 0.45. There is a generic 
underabundance of DM throughout the wino/higgsino dominated LSP region of 0200 ^ 0.15. 
On the other hand, 0:200 < 0.1 results in a dominant bino content of the LSP leading to a 
generic overabundance of DM. 

Thus we have found the optimised mixing parameters to be 

075 = 0.475 and 0200 = 0.12 (14) 

so as to get the most favourable DM relic density. There is admittedly some finer adjustment 
involved in these parameters so that we may be able to probe such large mixing zones of 
LSP. This is analogous to the delicate correlations between the mass parameters existing in 
the stau coannihilation or the resonant annihilation regions that provide with cosmologically 
favourable DM relic density in mSUGRA models, as pointed out in Section 1. 

Note however that once the mixing parameter is so adjusted, then the (1 + 75) and 
(1 + 200) models are each as predictive as the mSUGRA model. Note that in the first case 
we need roughly equal contributions from the singlet and the 75-plet superfields to achieve 
a favourable DM relic density, while in the second case it is achieved through a dominantly 
singlet contribution with only ~ 10% admixture from the 200-plet superfield. 

3 DM Relic Density for the (1 + 75) Model: 

We have computed the DM relic density for the (1 + 75) and (1 + 200) models in terms of the 
GUT scale mass parameters mo and mi/2 using SUSPECT [25] for the numerical evaluation 
of the two-loop RGEs. The DM relic density was evaluated using micrOMEGAs [26] and 
cross-checked with DARKSUSY [27]. 

Fig.l shows the DM relic density in the mo — mi/2 plane for the (1 + 75) model with the 
mixing parameter of eq.(14) and a moderate value of tan/? = 10. For simplicity we have 
set Ao = and positive /i. The upper disallowed region (in gray) corresponds to the lighter 
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Figure 1: Allowed and disallowed zones of the (1 + 75) model. The WMAP DM relic density 
satisfying regions are shown in red in the mo — mi/2 plane for tan/3 = 10, = 0, /i > 
corresponding to ays = 0.475 of Eq.l4. The upper disallowed region (in gray) corresponds to 
no EWSB (/i^ < 0). The lower gray region corresponds to stau becoming the LSP. Contours 
are drawn for Mi, fi, Zb and VLcDKih'^- 

chargino being smaller than the LEP limit [3] or no EWSB (/i^ < 0); and the lower one to a 
stau LSP. The allowed region is mapped by constant Mi and /i contours. At the intersection 
points Ml = yU, the LSP is an equi-mixture of bino and higgsino {Zb = c\ = 0.5). The bino 
fraction increases with decreasing mo, as indicated by the fixed Zb contours. The region 
satisfying the DM relic density value of the WMAP data [4] 

VtcDuh^ = Q-l^QQ ±0.0186 (3(t), (15) 

is indicated by the bands of red dots. The contours of fixed QcDAih^ = 0.05 and 0.2 are 
also shown on the two sides of these red bands. 

We see from Fig.l that almost half the parameter space corresponds to DM relic density 
being in the right ball-park of QcDMh"^ = 0.05 — 0.2. And even the precise WMAP value of 
Eq.(15) is satisfied by two fairly thick bands passing through the middle of the parameter 
space. The upper band corresponds to the DM being a roughly equal admixture of B — H, 
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Table 2: MSSM masses in GeV for a few sample parameter points for the 1 + 75 model. 



where their pair-annihilation occurs mainly via the gauge coupling of the higgsino compo- 
nent. With decreasing mo the bino component (Zb) increases, leading to the increase of 
DM relic density upto QcDAih"^ — 0.6. However further decrease of mo leads to resonant 
pair-annihilation via the pseudo-scalar higgs boson A, as 2Mi — > m^. This is the dominant 
dark matter annihilation mechanism for the lower red band, satisfying the WMAP relic den- 
sity of Eq. (15). The region below this band has an underabundance of DM relic density 
because of the resonant pair annihilation as m^ — 2Mi. This is the so-called funnel region. 
The lower edge of the funnel is marked by the red strip at the bottom right, adjacent to the 
lower boundary. We have found this strip to be dominated by resonant pair annihilation of 
DM via A, with only ~ 10% contribution from the stau coannihilation. The QcDAih'^ — 0.05 
contour (its lower branch) practically overlaps with the upper edge of this red strip. 

Table 2 lists the SUSY spectra for three representative points on each of the two main 
branches, satisfying DM relic density from WMAP(Eq.l5). Note that the SUSY spectra 
show an inverted mass hierarchy , where the lighter stop(ti) is much smaller than the other 
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squark masses. This is related to the large negative contribution from the terms associated 
with top- Yukawa coupling in the renormalisation group evolution of the third generation of 
squark masses [11]. This feature is more pronounced for the upper branch, represented by 
the first three columns due to larger values of mo- But it is quite significant for the lower 
branch as well. Note also an approximate degeneracy among the three lighter neutralinos 
(X123) ^he lighter chargino (Xi). This is again more prominent for the upper branch, 
which corresponds to a larger mixing between higgsinos and bino. Consequently, the coan- 
nihilation of Xi with xt Xi" make important contributions to the annihilation process. 
The dominant annihilation processes for the upper branch are 

XiXi, xlxl, Xixt ZZ, ff (both via s and t - channel processes) (16) 

which are driven by the gauge coupling of the higgsino component. On the other hand the 
lower branch is dominated by the resonant annihilation processes 

XiXi, XiX2^ XiXt ti, bb,Tf, tb, tp^ (via s - channel A, H, H+). (17) 

Since the gauge coupling of W/Z to the LSP pair goes like the square of its higgsino com- 
ponent, the annihilation process of Eq.l6 are strongly suppressed with the decrease of the 
higgsino component to ~ 10% for the lower branch. On the other hand the A coupling 
to the LSP pair goes like the product of its higgsino and gaugino components and hence 
remains significant for a higgsino component of ~ 10%. Even more importantly there is a 
large resonance enhancement for this region as 2m^o ttla- This also includes some amount 
of coannihilation of Xi and Xi the s-channel. 

Note finally that a higgsino component of ~ 10% in the LSP is still large compared to 
the mSUGRA model except a few regions like HB/FP zones. This is why one can get a 
resonant annihilation region here even for a moderate value of tan (3 [28] . 

4 Direct & Indirect Detection rates for the (1 + 75) 
Model: 

We have estimated the DM signals for the (1 + 75) model for direct and indirect detection 
experiments using the DARKSUSY [27]. The direct detection signal is based on the elastic 
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scattering of the DM on a heavy nucleus hke Germanium or Xenon. The main contribu- 
tion comes from the spin-independent interaction of the DM Xi with nucleon, which adds 
coherently in the nucleus. This is dominated by the Higgs exchange. Fig. 2 (a) shows the 
scatter plot of this cross-section against the DM mass. The discovery limits of present and 
proposed direct detection experiments are also shown for comparison. We have shown the 
limits from CDMS (Ge) 2005 [29], XENON-10 [30] and future SuperCDMS (Snolab) [31] 
and XENONIT [32] experiments in the figure. The upper and lower shaded (red) branches 
correspond respectively to the upper and lower WMAP relic density satisfying branches of 
Fig.l. Since the Higgs coupling to the neutralino DM is proportional to the product of its 
higgsino and gaugino components, one gets a larger cross-section for the upper branch of 
Fig. 2 (a), where DM corresponds to a roughly equal admixture of bino and higgsino. This 
comes within the detection limit of the proposed SuperCDMS experiment [31]. Both the 
branches will be covered by the proposed 1-Ton Xenon [32] experiment. Note that these 
cross-sections are larger than those of the mSUGRA model except the focus point region. 
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(a) (b) 
Figure 2: (a): Spin-independent neutralino-proton cross section for the (1 + 75) model 
vs. LSP mass for the parameters shown in Fig.l. The WMAP satisfied regions are shown 
in red. Various limits from the recent (DAMA, CDMS and XENON-10) and the future 
(SuperCDMS and XENON- IT) experiments are shown, (b): Spin-dependent neutralino- 
proton cross section for the (1 -|- 75) model vs. LSP mass for the parameters shown in Fig.l. 
The WMAP satisfied regions are shown in red. 



Fig. 2(b) shows the corresponding scatter plot of the spin-dependent cross-section, which 
is dominated by Z exchange. Again upper and lower shaded (red) branches correspond to 
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WMAP satisfying branches of Fig.l respectively. Since the Z couphng to the neutrahno DM 
is proportional to the square of its higgsino component, one gets a larger cross-section for 
the upper branch of Fig.l, corresponding to a larger higgsino component of DM. The spin- 
dependent cross-sections in this model are significantly larger than those of the mSUGRA 
model (except the focus point region). But still they are much below the detection limits of 
any present or proposed direct detection experiments [33]. 
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Figure 3: Neutralino annihilation induced muon flux from the Sun in km~'^yr~^ vs LSP mass 
for the (1 -|- 75) model for the parameters shown in Fig.l. The WMAP satisfied regions are 
shown in red. 



A very promising indirect detection experiment for DM is via the high energy neutrinos 
coming from their pair-annihilation in the solar core. Since the annihilation rate at equi- 
librium is balanced by the DM capture rate inside the Sun, the resulting neutrino signal is 
proportional to the Xi ~ P cross-section. This is dominated by the above-mentioned spin- 
dependent interaction via Z exchange. Fig. 3 shows the model prediction for the rate of muon 
signal events, resulting from these high energy neutrinos, in a km^ size neutrino telescope 
like IceCube [34]. Again upper and lower shaded (red) branches correspond to the upper 
and lower WMAP satisfying branches of Fig.l. One can even see a small red strip at the 
right end, corresponding to a similar one in Fig.l. One sees a very promising signal with 
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> 10 events/year at the IceCube. Since the Z couphng to the neutrahno DM goes hke the 
square of its higgsino component, the size of this signal is much larger than that of the bino 
dominated DM of the mSUGRA model. 

We have also computed the high energy 7-ray signal coming from the pair annihilation 
of the DM at the galactic core, assuming the standard NFW profile of DM distribution near 
the galactic core [35]. Fig. 4(a) shows the signal rate of line 7-rays coming from the DM pair 
annihilation, XX — > 77(7^) with ~ via VT-boson loop. Again this decay amplitude is 
proportional to the square of the higgsino component of the neutralino DM. Consequently 
the upper and lower shaded (red) branches correspond to the respective WMAP satisfying 
branches of Fig.l, while the small red strip at the right end correspond to the corresponding 
one of Fig.l. 

Fig. 4(b) shows the signal rate of continuum 7-rays coming from the tree level pair an- 
nihilation processes of Eq.l6 and Eq.l7. Production of neutral pions and their subsequent 
decays into photons is the most important mechanism for this signal [37]. While the 66 
channel is the most prominent channel for the continuum 7-rays signal, the WW and ZZ 
channels are also very significant. The two shaded (red) bands of Fig.4(b) correspond to the 
two WMAP relic density satisfying bands of Fig.l. The isolated red strip on the top right 
corresponds to that on the bottom right of Fig.l, which is dominated by the 66 channel. 
Finally Fig. 4(c) shows the continuum 7-ray spectrum for some representative points of this 
model, along with the discovery limit of the Fermi Gamma-ray Space Telescope (FGST) 
experiment or formerly known as GLAST [38,39]. One would need a boost^ factor of ~ 10^ 
to see this 7-ray signal in FGST/GLAST. Note however that one can raise (lower) the signal 
rate by a factor of ~ 10"^ by assuming a more spiked (flat) DM profile like that of Moore 
profile [41] (or a spherically symmetric isothermal core profile [42]). 

5 DM Relic Density for the (1 + 200) Model: 

Fig. 5 shows the DM relic density in the mo — vtixj^ plane for the (1 + 200) model with the 
mixing parameter of eq.(14). Again the upper and lower disallowed regions corresponds 

"'^Processes like adiabatic compression may cause a large enhancement of the DM density near the galactic 
center (see e.g. [36]). Larger DM annihilation rate may also result from existence of clumps within the halos 
of galaxies [40]. A boost factor essentially estimates a few such effects. 
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Figure 4: (a): Monochromatic 7-ray flux from DM pair anniliilation near tlie galactic center 
vs. the LSP mass shown for the NFW profile of DM halo distribution and an aperture size 
of AQ = 10~^ sr for the (1 + 75) Model. The parameters are same as in Fig.l. The WMAP 
satisfied regions are shown in red. (b): Same as (a) except that the vertical axis refers to a 
continuous 7-ray flux above a threshold energy of 1 GeV. (c): Photon spectra expected from 
two particular LSPs corresponding to mo = 270 GeV and mo = 1350 GeV for the (1 + 75) 

Model with mi = 400 GeV for an NFW profile. The other parameters are same as in Fig.l. 

2 

The upper set of curves corresponds to a boost [36] factor of 10^. The discovery limit of 
FGST/GLAST is shown blue straight line. 



no EWSB (yU^ < 0); and to a stau LSP respectively. The allowed region is mapped by 
constant Mi and M2 along with the constant /i contours. We see from these contours 
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that the LSP is an equimixture of bino and wino (B-W), along with a significant higgsino 
component in the mo > 1 TeV region. In fact the intersection points of these contours 
(Ml, M2, /i = 500 and 800 GeV) mark the region of equally mixed B-W-H LSP. The 
combined gaugino component of the LSP is indicated by the fixed Zg{= c\ + c^) contours. 
This, when subtracted from unity would be a measure of the higgsino component. The 
region satisfying the DM relic density value of Eq.l5 from WMAP is indicated by the bands 
of red dots. We see two thick bands in the middle of the parameter space. Note that these 
DM relic density satisfying bands correspond to LSP mass > 500 GeV. The contours of fixed 
^CDAih"^ = 0.05 and 0.15 are also shown on the two sides of these red bands. The bulk of the 
parameter space outside these red bands corresponds to underabundance of DM relic density. 
In particular the corridor cutting through the two bands correspond to underabundance due 
to rapid resonant annihilation of DM pair via A-boson (i.e. the funnel region). For a clearer 
understanding of this model we list the SUSY spectra for two representative points from each 
of the two DM relic density satisfying branches in Table 3. This model requires a rather heavy 
SUSY spectrum for satisfying the DM relic density constraint. Note that the SUSY spectra 
of this model show an inverted mass hierarchy like the previous model, where the lighter stop 
(ti) is significantly lighter than the other squarks. Note also the close degeneracy among the 
lighter neutralino and chargino states (x? 2 ^"^^ xt) large bino- wino mixing. Therefore 

the coannihilation of Xi with X2 xf make important contributions to the annihilation 
process as in the previous case. Both the branches are dominated by the annihilation process 
of Eq.l6, driven by the gauge coupling of the wino (along with that of higgsino at large mo). 
But there is a 30 — 40% contribution coming from the resonant annihilation processes of 
Eq.l7 on the strips adjacent to the 'funnel' corridor, where 2m^o rriA. Note that the 
higgsino component of the LSP in this region is > 10%, which is large enough for resonant 
annihilation via A even for a moderate value of tan/? (= 10). 



6 Direct & Indirect Detection rates for the (1 + 200) 
Model: 

Fig. 6(a) shows the spin-independent scattering cross-section of the DM on nucleon, which 
gives the direct detection signal. The shaded (red) region covers the two WMAP relic 
density satisfying branches of Fig. 5. Most of this region corresponds to the upper branch 
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m,,,(GeV) 



Figure 5: Allowed and disallowed zones of the (1 + 200) model. The WMAP DM relic density 
satisfying regions are shown in red in the mo — mi/2 plane for tan/3 = 10, = 0, /i > 
corresponding to 0200 = 0.12 of Eq.l4. The upper disallowed region (in gray) corresponds to 
no EWSB (/i^ < 0). The lower gray region corresponds to stau becoming the LSP. Contours 
are drawn for Mi, M2, fi, Zg and QcDAih"^- 

of Fig.5, except for the lower edge, which corresponds to the lower branch. There is a thin 
line separating the two regions, which is barely visible in this figure. Note that the spin- 
independent scattering is dominated by the Higgs exchange, whose coupling to the neutralino 
DM is proportional to the product of its higgsino and gaugino components. Thanks to the 
presence of 10 — 20% of higgsino component in this model (Fig.5), the spin-independent 
cross-section is large enough to be detectable at the proposed superCDMS experiment [31]. 
Fig. 6(b) shows the model prediction for the spin-dependent cross-section of the DM on 
nucleon. This is dominated by the Z exchange whose coupling to the neutralino DM is 
proportional to the square of its higgsino component. The shaded (red) region covers the 
two WMAP relic density satisfying branches of Fig.5. Again most of this region corresponds 
to the upper branch, with only the lower edge corresponding to the lower branch. As in the 
case of the previous model, the spin-dependent cross-sections here is too small to be seen at 
direct detection experiments [33]. 
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parameter 


A 




B 


C 


D 




725 




900 


725 


900 


mo 


1450 




1357 


590 


950 




792 




983 


846 


1009 


Ml 


652 




814 


646 


811 




633 




798 


633 


797 


X2 


657 




818 


650 


815 




794 




985 


848 


1011 




822 




1009 


869 


1032 




643,81 


8 


807,1005 


641,865 


806,1028 




1700 




2045 


1637 


2017 




1460 




1649 


1216 


1540 




1813,1 


801 


2013,2001 


1478,1452 


1860,1843 




1535- 


1555 


1505-1530 


800-830 


1160-1190 




2000- 


2050 


2170-2240 


1520-1600 


1950 -2040 


"^a(^ mH+,mH) 


1726 




1799 


1174 


1546 



Table 3: MSSM masses in GeV for a few sample parameter points for the 1 + 200 model. 

Similar to what was discussed in Section 4 a larger Higgsino content of the LSP for 1 + 200 
scenario makes indirect detection of high energy neutrinos interesting. Fig. 7 shows the rate 
of the resulting muon events for the IceCube experiment. The upper and lower shaded (red) 
branches correspond to the respective WMAP satisfying branches of Fig. 5, with a clearly 
visible separation between them. As in the previous model, one expects a detectable signal 
rate > 10 events per year at the IceCube experiment. 

Fig. 8 shows the rates of the line and the contunuum 7-ray signal coming from DM pair 
annihilation at the galactic core assuming the NFW profile of DM distribution [35]. Fig. 8 (a) 
shows the signal rate of the line 7-rays coming from DM pair annihilation, xx 77(7^) 
via W boson loop. Since the wino has a large Isospin gauge coupling, the resulting 7 ray 
signal here is at least an order of magnitude larger than the previous model. The shaded 
(red) band corresponding to the WMAP satisfying region of this model, predicts a line 7-ray 
flux of ~ 10"^^ cm~^ s~^. Several atmospheric Cerenkov experiments have reported 7-ray 
events [43-46] at higher rate than this, but with continuous power-law energy spectrum. 
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(a) (b) 
Figure 6: (a): Spin-independent neutralino-proton cross section for the (1 + 200) model 
vs. LSP mass for the parameters shown in Fig.5. The WMAP satisfied regions are shown 
in red. Various hmits from the recent (DAMA, CDMS and XENON-10) and the future 
(SuperCDMS and XENON-IT) experiments are shown, (b): Spin-dependent neutrahno- 
proton cross section for the (1 -|- 200) model vs. LSP mass for the parameters shown in 
Fig.5. The WMAP satisfied regions are shown in red. 

typical of supernova remnants. Therefore it will be hard to separate the DM signal from 
this background with the present energy and angular resolution of the ACT experiments. 

Fig. 8(b) shows the signal rate of the continuum 7-rays, coming from the tree- level pair 
annihilation processes of Eqs.l6 and 17, with red dots corresponding to the WMAP sat- 
isfying regions. Since the most copious source of these 7 rays is the 6-quark jet from Eq. 
17, the upper red dots correspond to the region of Fig.5 with large resonant annihilation 
contribution. Finally Fig. 8(c) shows the continuum 7-ray spectrum for some representative 
points of this model, along with the discovery limit of the FGST/GLAST experiment [47]. It 
shows that a boost factor of ~ 10^ to see this 7-ray signal in FGST/GLAST for the assumed 
NEW profile of DM distribution. 

7 LHC Signatures of the (1 + 75) and (1 + 200) Models: 

A quantitative investigation of the LHC signature of these two nonuniversal gaugino mass 
models is beyond the scope of this work. We shall only remark on some prominent features, 
which are quite evident from the SUSY spectra of Table 2 and 3. Table 2 clearly shows an 
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Figure 7: Neutralino annihilation induced muon flux from the Sun in km~'^yr~^ vs LSP mass 
for the (1 + 200) model for the parameters shown in Fig.5. The WMAP satisfied regions are 
shown in red. 

inverted mass hierarchy for the (1 + 75) model where the lighter stop is about 30% lighter 
than the 1st and 2nd generation squarks. Moreover it is lighter than the gluinos over most 
of the parameter space. So one expects a SUSY signal from the channels 

pp ^ gg ^ ttiiii, pp^tiii (18) 

inclusive of their antiparticles. In view of the singlet dominance of ti, one expects it to decay 
via mainly gauge (Yukawa) coupling to B{H) into 

h ^ 2,3, bxt. (19) 

Since all these neutralino and chargino states are nearly degenerate, we expect their p^ to 
be largely carried by the LSP(xi), giving relatively large missing-p^ {f^r) events along with 
soft jets (leptons). Note however that one expects hard leptons and 6-jets from the decay 
of the top quarks in (18) and (19), resulting in hard and isolated multilepton+missing pj- 
signal accompanied by 4(2) 6-jets for the pair production of gluon (stop). 

The SUSY spectra of Table 3 also shows an inverted hierarchy for the (1+200) model with 
the ti being much lighter than the other squarks. Moreover it is lighter than the gluino over 
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Figure 8: (a): Monochromatic 7-ray flux from DM pair anniliilation near tlie galactic center 
vs. the LSP mass shown for the NFW profile of DM halo distribution and an aperture size 
of An = 10~3 sr for the (1 + 200) Model. The parameters are same as in Fig. 5. The WMAP 
satisfied regions are shown in red. (b): Same as (a) except that the vertical axis refers to a 
continuous 7-ray flux above a threshold energy of 1 GeV. (c): Photon spectra expected from 
two particular LSPs corresponding to mo = 750 GeV and niQ = 1225 GeV for the (1 + 200) 
Model with mi = 800 GeV for an NFW profile. The other parameters are same as in Fig.5. 
The upper set of curves corresponds to a boost [36] factor of 10^. The discovery limit of 
FGST/GLAST is shown blue straight line. 



essentially the entire parameter space. Because of the rather heavy mass range of squarks 
and gluino in this model we expect a very important contribution to the LHC SUSY signal 
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to come from the pair production of ti. Note that in this case the dominant decay channels 
are 

^1 ^ ^X?,2,3,4, bxt, (20) 

since the higgsino dominated states are the heavier chargino and neutrahno states {xt , Xsa)- 
However these states are also quite close to the LSP Xi this case. Therefore one expects 
signal characteristics as the earlier model. It should be mentioned here that similar char- 
acteristics of the LHC signal was noted for the focus point region of the mSUGRA model 
in [11], and has since been studied in detail in [12]. One expects similar features for these 
nonuniversal gaugino mass models of the mixed neutralino LSP as well. Only in these models 
one can have relatively light scalar mass mo, so that the lighter ti is seen to be even lighter 
than the gluino, unlike the focus point case. 

8 Conclusion 

We have investigated the dark matter properties of non-universal gaugino mass scenarios 
where the gauge kinetic energy function is a mixture of chiral superfields that transforms 
as singlet and the 75-dimensional or the singlet and the 200 dimensional representation 
of SU(5). The mixing of the representaions are chosen so as to probe the dark matter 
annihilation properties of LSPs that are either a strong mixture of bino and Higgsinos in the 
1+75 model or a mixture of bino-wino or even a mixture of bino-wino-higgsinos in the 1+200 
model. For each of the scenarios we have analyzed the dark matter relic density in relation 
to the WMAP data and we have obtained broad regions of parameter space in the mo — mi 

2 

plane that satisfy the WMAP data. We have identified the annihilation and coannihilation 
channels. We have further computed the direct detection rates of neutralino DM for the 
spin-independent and spin- dependent Xi ^ P cross sections. The spin-independent cross 
section ranges would be successfully probed in the near future via XENON-IT. We have 
also computed the neutralino-annihilation induced muon flux from the Sun that originates 
from the high energy neutrinos produced via pair-annihilation of LSPs in the solar core. 
We have found promising signals that will easily be probed in the IceCube experiment. We 
have also obtained interesting photon signal rates (both monochromatic and continuum) for 
the indirect detection of LSPs via 7-rays. FGST/GLAST will be able to probe a significant 
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region of parameter space with a possible boost factor. Finally, we have described a few 

distinctive features of the expected LHC signals of these models. 
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